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The algorithm is described that enables one to perform an explicit summation of all the 
(n 2 / ln(Q 2 /A 2 )) JV corrections to a s (Q 2 ) that appear owing to the analytic continuation from 
spacelike to timelike region of momentum transfer. 



I. INTRODUCTION 



Perturbative QCD is intensively applied now {1} to various processes involving large momentum transfers, both 
in spacelike (q 2 = —Q 2 < 0) and timelike (q 2 > 0) regions (for a review see (||). However, the coupling 
constant g(n) (i.e., the expansion parameter) is defined usually with the reference to some Euclidean (spacelike) 
configuration of momenta of scale \i. For spacelike q this produces no special complications. One simply uses the 
rcnormalization group to sum up the logarithmic corrections (<? 2 (m) ln(Q 2 / /j, 2 )) n that appear in higher orders and 
arrives at the expansion in the effective coupling constant a s {Q 2 ) which in the lowest approximation is given by 
the famous asymptotic freedom formula E| 

2 47T 

as{Q ) = (11 -2iV / /3)ln(Q 2 /A 2 ) ' (1) 

where A is the "fundamental" scale of QCD. In general, the A-parametrization of a s (Q 2 ) is a series expansion in 
1/L (where L = ln(Q 2 /A 2 ) ), and the definition of A is fixed only if the 0(l/L 2 )-term is added to Eq.(l) j|. 

For timelike q there appear, however, i7r-factors (ln(Q 2 //i 2 ) — > ln(Q 2 / /i 2 ) ±wr), and it is not clear a priori what is 
the effective expansion parameter in this region. This problem was discussed recently {1} in a very suggestive paper 
by Pennington and Ross ||. These authors analyzed the ratio R(q 2 ) = <r(e + e~ — > hadrons)/cr(e + e~ — » fi + fi~) for 
which the analytic continuation from the spacelike to timelike region is well-defined and investigated which of the 
three ansatze (a s (q 2 ), \a s (—q 2 )\ and Kea s (—q 2 ) ) better absorbs the (7r 2 /L 2 ) JV -corrections[] in the timelike region 
q 2 > 0. Their conclusion was that \a s (—q 2 )\ is better than a s (q 2 ) and Rea s (— q 2 ). Nevertheless, it is easy to 
demonstrate by a straightforward calculation that \a s (—q 2 )\ cannot absorb all the (tt 2 / L 2 ) N -terms associated with 
the analytic continuation of the ln((5 2 //i 2 )-factors. Our main goal in the present letter is to show that by using 
the A-parametrization for a s (Q 2 ) in the spacelike region it is possible to construct for R(q 2 ) in the timelike region 
the expansion in which all the (tt 2 / L 2 ) N -terms are summed explicitly. 

II. A-PARAMETRIZATION IN SPACELIKE REGION 



The starting point for the A-parametrization is the Gell-Mann-Low equation taken as a series expansion in 
G = a s /4n: 

L ee ln(Q 2 /A 2 ) = J_ + JilnG + A + Q + {G 2 ) , (2) 

where bk are /3-function coefficients: 

&o = 11 - 2N f /3 §, bx = 102 - 38^//3 §, bf s = 2857/2 - 5033AT//18 + 325^/54 §. 

The parameter A in Eq.(2) is due to the lower boundary of the GML integral (8]||. By a particular choice of A 
one fixes the definition of A: A = A(A) Eq.(2) is solved by iterations and the result is reexpanded in 1/L: 

where 

Li = % HboL) - A . (4) 

The expansion (3) is useful, of course, only if it converges rapidly enough. In fact, the convergence of the l/L 
series depends (i) on the value of L we are interested in and (ii) on the choice of A. 

We emphasize that the most important for perturbative QCD is the region L > 3, since L = 3 corresponds to 
a s ~ 0.5, and the reliability of perturbation theory for larger a s is questionable. Hence, in a realistic situation the 
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(3) 



*Odd powers of {in / L) cancel because R is real. 

^Of course, A depends also on the renormalization scheme chosen. 



2 



naive expansion parameter 1/L is smaller than (but usually close to) one third. Of course, 1/3 is not very small, 
so one must check the coefficients of the 1/L expansion more carefully. First, there is a A-convention- independent 
term (62^0 — ^/(b^L 2 ) which reduces for Nf — 3 to roughly 0.25/L 2 and gives, therefore, less than 3%-correction 
to the simplest formula (1). There are also A-dependent terms like Li/L, Li/L 2 and one should choose A so as to 
minimize the upper value of the ratio L\/ L in the L-region of interest. 
If one takes, e.g., 

A = A opt = (6 1 /^)l n (46 ) 

then L\ = (&1/&0) ln(L/4) and the ratio L\/L is smaller than 7% in the whole region L > 3. Another choice |l(]] 
is to take 

& = A(Q 2 ) = (b 1 /b 2 )]n(b L ), 

where Lo = Ih(Qq/A 2 ) and Q\ lies somewhere in the middle of the Q 2 -region analyzed. In this case 
L\ = {bi/b 2 )) Ih(L/Lq), i.e., L\/L is zero for Q 2 ~ Qq and smaller than 7% for all Q 2 in the region where L > 3. 
An important observation is that both the choices minimize the corrections not only in Eq. (3) but also in the 
GML equation (2). 

Really, for small G the only dangerous term in Eq. (2) is In G, hence, the best thing to do is to compensate it by 
taking A = — (&i/&q) In G, where G is a s (Q 2 )/47r averaged (in some sense) over the relevant Q 2 -region. After this 
has been done, one may safely solve Eq. (2) by iterations and perform the 1/L-expansion. For a proper choice of 
A Eq. (3) has 1% accuracy for L > 3, and, moreover, the total correction to the simplest formula (1) is less than 
10%. However, accepting the most popular prescription 

A pop = (&i/&g)ln(&o) = A(Q 2 = eA 2 ) 

(the only motivation for A pop being the "aesthetic" criterion that L x should have the shortest form 
L\ = (&i/6q) ln(-L) ) one minimizes L\/L in the region Q 2 ~ 3A 2 nobody is really interested in. Moreover, in 
the important region L ~ 3 one has L^ op /L ~ 1/3, and the convergence of the 1/L series is very poor in this case. 

Thus, the A-parametrization (Eq.(3)) gives a rather compact and sufficiently precise expression for the effective 
coupling constant in the spacelike region provided a proper choice of the A-parameter has been made. 



III. A-PARAMETRIZATION AND R(e+e -> hadrons, s) 



The standard procedure (see, e.g., |jTl| and references therein) is to calculate the derivative D(Q 2 ) — Q 2 dt/dQ 2 
of the vacuum polarization t(Q 2 ) related to R by 

R(s) = ^(t(-s + ie)-t(-s-ie)) . (5) 

In perturbative QCD D{Q 2 ) is given by the a s (Q 2 )-expansion: 



Only di is known now |ll| , [12| {2}, its value depending on the renormalization scheme chosen. Using Eq. (5) and 
the definition of D, one can relate R(s) (or, more precisely, its perturbative QCD version R^ CD (s) ) directly to 
D(Q 2 ) 

R QCD (s) = ^-. D(a)™. (7) 

Integration in Eq.(7) goes below the real axis from — s — it to zero {3} and then above the real axis to — s + it. 
In a shorthand notation D => R = $[£>]. In some important cases the integral can be calculated explicitly {4} : 



1^1 (8) 
- ^-arctan (7 r/L s ) = 1 



arctan(7r/L s ) = —!--— + ... (9) 
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\n{L a /L ) ln(Vg + 7r 2 /L ) - (L s /tt) arctan(7r/£,) + 1 _ L s /L L_* 



T 2 



5 7T^ 
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i» * <_1) "(^ij! (sl:) zrn? = i» { 1- ii ( e ■ + " •} (12) 

where L s = ln(s/A 2 ), L CT = ln(er/A 2 ) and Lo is a constant depending on the A-choice. 

Using the A-parametrization for a s (a) and incorporating Eqs.(8)-(12) (as well as their generalizations for 
In 2 L/L 3 , \nL/L 3 etc.) produces the expansion for R9 CD {s) 



Jl I t2\N 



q I fe=l J 



(13) 



in which all the (n / L ) -terms are summed up explicitly. 

IV. QUEST FOR THE BEST EXPANSION PARAMETER 



Note that the expansion (13) is not an expansion in powers of some particular parameter since the application of 
the ^-operation normally violates nonlinear relations: $[1/L 2 ] ^ ($[1/L]) 2 , etc. A priori, there are no grounds to 
believe that a power expansion is better than any other (say, Fourier). In fact, the expansion (13) converges better 
than the generating expansion (6) for D(a) because, as it follows from Eqs. (9)-(12), $[0!^] is always smaller than 
. Moreover, (3>[a^ r+1 ]) 1 '( N+1 ' < (^[a^]) 1 ^ , i.e., the effective expansion parameter decreases in higher orders. 
Thus, if one succeeded in obtaining a good expansion for D(a) (with all cLm being small numbers), then the 
resulting ^af^] expansion for R9 CD {s) is even better, and the best thing to do is to leave it as it is. 

However, if one insists that the result for R9 CD {s) should have a form of a power expansion, then the best 
expansion parameter is evidently <I>[a s /7r] because the largest nontrivial (i.e., 0(a s /7r) ) term of the expansion is 
reproduced in the exact form and only higher terms are spoiled. The analogue of the simplest A-parametrization 
for a s (Q 2 ) (Eq.(l) ) is then 

a s (g 2 ) = A arctam (_!L_) . (14 ) 

Using Eqs. (8)-(13) it is easy to realize that a s (q 2 ) is really a bad expansion parameter, because if one reexpands 
a s (q 2 ) in a s (q 2 ), there appear terms with large coefficients 




If one reexpands a s (q 2 ) in Re a s (— q 2 ) then the corresponding coefficient is even 2 times larger, whereas if a s (q 2 ) 
is reexpanded in \a s (— q 2 )\, the coefficient is 2 times smaller. This observation is in full agreement with the result 
of ref. H quoted in the Introduction. 



V. CONCLUDING REMARKS 



It should be noted that the change of the expansion parameter as given by Eq. (15) affects only the (a s /7r) 3 
coefficient of the f^^-expansion which has not been calculated yet {2}. So, within the present-day accuracy, 
all expansions for RP CD have the same coefficients. It is worth emphasizing, nevertheless, that the tt 2 /L 2 terms 
produce for a s > 0.3 more than 20% correction to a s , i.e., they are more important (for an optimal choice of the 
A-parameter) than the 2-loop corrections in Eq. (3). 

To conclude, we have described the construction of an optimized (i.e. rapidly convergent) A-parametrization for 
the effective QCD coupling constant in the spacelike region, and then we used it to obtain the fastest convergent 
expansion for the timelike quantity R^ CD (s). The technique outlined in the present paper can be applied also 
to other R < ^ CD -\ike quantities. Such quantities do appear, e.g., in the QCD sum rule approach |ll| in which the 
analysis of hadronic properties is based on the study of vacuum correlators of various currents. They appear also 
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in an alternative approach [|L4| based on the finite-energy sum rules |15[. It should be stressed that in the latter 
approach the RP CD -\ike quantities enter into the basic integral relation, and the analysis is most conveniently 
performed if one has a simple analytic expression similar to that described above. 
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NOTES ADDED 



{1} This paper was submitted in 1982 to Physics Letters B, but not accepted because the referee was not 
convinced that it needs a rapid publication. I was recommended to write a longer version and submit it to a 
regular journal. Unwisely, I did not do that. Still, though the paper existed in the preprint form only, it was 
known to experts in multiloop calculations. In particular, Eq.(15) was incorporated into the 4-loop calculation of 
RP (s) [16,17]. Later, some of my results were used or rediscovered in several publications including very recent 
ones (see, e.g., [18-22]). In 1996, the paper was reprinted in JINR Rapid Communications [23], but the subsequent 
experience convinced me that the only way to make the paper accessible to interested readers (and hopefully get 
credit for its results) is to submit it to an e-print archive. The original version is reproduced above without any 
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changes. I only added references to short notes presented below. They contain updating and clarifying remarks. 

{2} The coefficient d^ IS was calculated in refs. [16,17]. Starting from this level, the coefficients r k of the 
a s (s)-expansion for BP CD (s) 

^%) = E^{i + ^ + .- 2 (^) 2 + -- 3 (^) i + ... 

differ from dk- according to our Eq.(15) r 3 = d 3 — (7rfe ) 2 /48. 

{3} This is an incorrect description of the actual straightforward procedure which I used to get results displayed 
in Eqs. (8)-(12). The central idea of the paper is to represent D(a) as a sum of terms for each of which the integral 
(7) can be calculated as an explicit analytic expression like In In a and then simply take the difference of these 
integrals at — s + it and — s — it using the standard prescription that the cut of lnz is on the negative real axis. 
This is precisely what one should do to analytically continue t(Q 2 ) from the deep spacelike region. However, for 
terms in -D(cx) containing poles at a = A 2 this prescription is equivalent to integration from —s — ie below the real 
axis to some real point <tq > A 2 (rather than zero) and then above the real axis to — s + ie. 

{4} In 1982, the applications of perturbative QCD in the < s < A 2 region were not treated as reliable, so 
s > A 2 (i.e. L s > 0) is implied in Eqs. (8)-(12). Furthermore, the results are presented in a form most suitable 
for the 1/L S expansion. To reconstruct the original expressions valid both for positive and negative L s one should 
change 

arctan(7r/L s ) — * tt/2 — arctan(L s /7r) 
in Eqs. (8)-(12), but this form is not very illuminative for large L s . 
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